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The  design  of  terminal  guidance  law  with  impact  angle  constraint  is  required  for 
air-to-ground  guided  weapons  to  increase  their  warhead  effect.  The  variable  structure 
guidance  law  that  consists  of  diving  plane  guidance  and  turning  plane  guidance  equations 
with  impact  angle  constraint  is  derived,  and  the  saturation  function  is  introduced  into  the 
design  of  reaching  law  control  to  weaken  the  chattering  of  the  guidance  system.  The  influ¬ 
ence  of  four  guidance  parameters  (i.e.,  reaching  law  factor,  switching  item  gain,  angle  error 
item  factor,  and  boundary  layer  thickness)  on  guidance  performance  is  studied  and  three 
typical  constraints  (i.e.,  heating  rate,  normal  load  factor,  and  dynamic  pressure)  are 
analyzed.  An  optimization  model  is  established  for  this  problem  and  the  feasibility  of 
on-line  optimization  on  guidance  law  parameters  by  the  Sequential  Quadratic  Program¬ 
ming  (SQP)  algorithm  is  discussed  as  well.  Simulation  results  show  that  the  on-line  optimi¬ 
zation  of  the  derived  guidance  law  not  only  satisfies  specified  constraints,  but  also 
minimizes  the  fuel  cost  during  the  flying  course.  Moreover,  the  optimization  process  can 
be  completed  in  a  few  seconds  so  that  it  is  suitable  for  on-board  applications. 

©  2013  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

A  guidance  system  is  a  device  or  a  group  of  devices  used  to  navigate  a  ship,  aircraft,  missile,  rocket,  satellite,  or  other 
crafts.  It  generates  guidance  laws  which  take  input  from  the  navigation  system  and  use  targeting  information  to  send  signals 
to  the  flight  control  system  that  will  allow  the  vehicle  to  reach  its  destination. 

The  main  objective  of  guidance  systems  for  ground  striking  weapons  is  to  generate  suitable  commands  that  produce  zero 
terminal  miss-distances.  In  some  cases,  however,  many  lands  of  guided  weapons  are  expected  not  only  to  get  a  minimum 
miss-distance  but  also  to  achieve  a  desired  impact  angle,  so  that  the  warhead  of  the  weapons  can  acquire  better  kill  effect. 
A  typical  example  is  the  design  of  a  ‘homing’  system  for  certain  re-entry  vehicles.  The  terminal  attitude  angle  should  be  with¬ 
in  a  permissible  range  in  addition  to  the  miss-distance  requirement.  Another  example  is  the  guidance  of  air-to-ground  weap¬ 
ons  with  directive  warheads.  During  the  final  attack,  the  impact  angle  should  have  a  preferred  value  to  achieve  a  proper 
penetration  and  ensure  high  killing  probability  for  some  missiles  and  torpedoes,  which  attack  targets  with  heavy  armors. 
For  the  past  few  decades,  a  variety  of  methods  for  guidance  laws  with  terminal  impact  angle  constraints  have  been  exten¬ 
sively  studied.  Based  on  different  theories,  the  existing  guidance  laws  are  categorized  into  optimal  guidance  laws,  proportional 
navigation  guidance  laws,  variable  structure  guidance  laws  and  other  guidance  laws. 
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Optimal  Control  Theory  (OCT)  plays  a  significant  role  in  the  study  of  impact-angle-constrained  terminal  guidance.  Kim  and 
Grider  [1]  first  derived  the  optimal  and  suboptimal  guidance  laws  with  constraints  on  impact  angle,  which  seems  to  be  a 
pioneering  research  in  this  area.  Based  on  this,  many  studies  are  started.  A  planar  engagement  problem  with  terminal  impact 
angle  constraint  was  formulated  as  a  numerical  optimization  problem  in  Song  et  al.  [2],  wherein  the  angle  constraint  was 
treated  as  a  penalty.  Ryoo  et  al.  [3]  investigated  the  closed-form  solutions  of  optimal  guidance  laws  and  expressed  the  entire 
missile  states  using  polynomial  functions  in  terms  of  time-to-go.  At  the  same  time,  a  generalized  formulation  of  energy  min¬ 
imization  optimal  guidance  law  for  constant  speed  missiles  was  proposed  in  Lee  et  al.  [4]  to  control  impact  angle,  terminal 
acceleration  and  miss-distance.  Ryoo  et  al.  [5]  presented  the  optimal  guidance  laws  for  arbitrary  missile  dynamics  and  de¬ 
rived  the  specific  forms  of  the  state  feedback  for  lag-free  and  first-order  lag  systems.  Based  on  this,  a  new  optimal  guidance 
law  for  constant  speed  missiles  was  proposed  in  Ryoo  et  al.  [6]  by  including  a  time- varying  weighing  function  in  the  cost  of 
control  energy  to  improve  the  performance  of  the  guidance  system  against  external  disturbances  or  sensor  measurement 
noises.  For  further  the  research  of  Lee  et  al.  [4],  a  suboptimal  guidance  with  the  terminal  constraints  on  impact  angle,  accel¬ 
eration,  and  time-derivative  of  acceleration  was  presented  in  Lee  et  al.  [7],  wherein  the  control  input  consists  of  an  additional 
feed-forward  command  and  an  optimal  feedback  command.  In  Subchan  [8],  the  problem  of  air-to-ground  missile  guidance 
with  impact  angle  constraint  was  reinterpreted  using  OCT  and  transformed  into  a  minimum  integrated  altitude  problem, 
which  was  solved  by  an  indirect  method.  Different  from  this,  a  three-dimensional  guidance  law  with  impact  angle  constraint 
was  proposed  by  Ma  et  al.  [9]  using  Lyapunov  stability  theory.  Kim  et  al.  [10]  further  dealt  with  the  guidance  and  control 
system  to  impact  a  target  with  a  desired  impact  angle  for  precision  guided  bombs  and  derived  the  control  based  on  the  solu¬ 
tion  of  the  linear  quadratic  optimal  control  problem.  Then  a  nonlinear  suboptimal  guidance  law  was  synthesized  in  Oza  and 
Padhi  [11,12]  using  recently  developed  Model  Predictive  Static  Programming  (MPSP)  for  successful  interception  of  ground 
targets.  The  main  feature  of  this  guidance  law  is  that  it  accurately  satisfies  terminal  impact  angle  constraints  in  both  azimuth 
and  elevation,  simultaneously.  As  a  new  trial,  Xing  and  Chen  [13]  considered  the  constraints  on  both  terminal  fight  path  an¬ 
gle  and  terminal  angle  of  attack,  and  derived  the  guidance  law  by  solving  the  linear  quadratic  terminal  control  problem. 
More  recently,  the  analytic  solutions  of  the  generalized  and  optimal  impact-angle-control  guidance  laws  for  a  first-order 
lag  system  were  investigated  in  Lee  et  al.  [14,15],  respectively,  to  examine  the  effects  of  system  lag  on  performance  of 
the  guidance  laws.  Furthermore,  the  work  presented  in  [16]  illustrated  the  optimality  of  linear  time-varying  guidance  laws 
for  controlling  impact  angles  as  well  as  terminal  miss-distances  using  the  inverse  problem  of  OCT. 

Proportional  Navigation  is  another  approach  adopted  by  many  researchers  in  the  design  of  terminal  impact-angle-con- 
strained  guidance  laws  except  for  OCT.  It  is  well  known  that  the  standard  Proportional  Navigation  Guidance  (PNG)  law  can¬ 
not  achieve  an  acceptable  performance  in  terms  of  the  impact  angle.  It  is  necessary  to  find  new  ways  to  take  the  terminal 
angular  constraint  into  consideration  in  the  deviation  of  the  terminal  PNG  laws.  The  guidance  method  discussed  in  Kim 
et  al.  [17]  added  a  time-varying  bias  term  in  the  conventional  PNG  to  achieve  the  final  impact  angle  condition,  which  made 
assumptions  of  constant  velocity  and  small  error  angles.  Then  a  variable  structure  proportional  navigation  law  for  the  pas¬ 
sive  homing  missile  was  designed  by  Song  and  Zhang  [18]  to  hit  the  target  with  a  minimum  miss-distance  and  desired  im¬ 
pact  attitude  angle  within  the  required  overload.  Besides,  a  kind  of  Biased  Proportional  Navigation  Guidance  (BPNG)  to  get  an 
arbitrary  impact  angle  against  stationary  target  was  proposed  in  Jeong  et  al.  [19]  and  the  closed  form  solution  of  this  guid¬ 
ance  was  derived.  Furthermore,  an  adaptive  proportional  guidance  law  with  terminal  angular  constraint  was  derived  in  Hu 
and  Cai  [20]  for  the  ground  stationary  target  and  an  adaptive  guidance  parameter  updating  method  was  presented.  In  addi¬ 
tion,  proportional  navigation  based  guidance  laws  for  impact  angle  constrained  interception  of  moving  and  stationary  tar¬ 
gets  were  studied  in  Ratnoo  and  Ghose  [21-23].  Similar  as  these  research,  Erer  and  Merttopcuoglu  [24,25]  devised  a  method 
to  control  the  impact  angle  in  an  indirect  manner  by  employing  biased  Pure  Proportional  Navigation  (PPN)  as  the  guidance 
law,  and  then  the  impact  angle  value  against  a  stationary  target  became  a  function  of  both  initial  angular  conditions  and 
integral  of  the  bias  value.  Moreover,  a  guidance  law  that  is  capable  of  achieving  a  wide  range  of  impact  angles  was  designed 
in  Akhil  and  Ghose  [26]  based  on  the  BPNG,  which  used  small  angle  approximations  to  derive  the  bias  term  for  impact  angle 
requirement. 

In  the  past  few  decades,  Variable  Structure  Control  (VSC)  and  associated  Sliding  Mode  Control  (SMC)  have  been  developed 
and  applied  enormously.  Considering  their  good  performance  on  controlling  nonlinear  plants,  the  VSC  and  SMC  technologies 
have  been  introduced  into  the  design  of  guidance  laws  to  improve  their  robustness.  In  early  1990,  the  VSC  theory  was  ap¬ 
plied  to  guidance  law  design  of  air-to-air  missile  by  Brierley  and  Longchamp  [27],  following  which  numerous  guidance  laws 
emerged  based  on  different  background.  A  homing  guidance  with  terminal  angular  constraint  was  stated  in  Kim  et  al.  [28]  by 
employing  a  Lyapunov-like  function  with  the  SMC  methodology,  but  the  problem  is  how  to  choose  suitable  guidance  param¬ 
eters.  Then  a  new  guidance  method  called  optimal  sliding  mode  guidance  law,  which  integrates  optimal  guidance  with  slid¬ 
ing  mode  guidance  was  derived  by  Zhou  et  al.  [29]  to  achieve  robustness  against  target  maneuver.  And  then  an  adaptive 
reaching  law  of  the  sliding  mode  was  presented  and  used  to  derive  an  adaptive  sliding  mode  guidance  law  by  Zhou  et  al. 
[30],  Based  on  these  studies,  a  scheme  of  integrated  guidance/autopilot  design  for  missiles  with  impact  angle  constraints 
was  studies  in  Xu  et  al.  [31  ]  based  on  the  back-stepping  idea  and  VSC  theory.  And  a  new  three-dimensional  variable  structure 
guidance  law  was  discussed  in  Sun  and  Zheng  [32]  based  on  adaptive  model-following  control.  For  further  understanding 
and  the  research  of  guidance  laws,  Hu  et  al.  [33]  derived  a  kind  of  variable  structure  guidance  law  with  terminal  angular 
constraint  and  discussed  the  feasibility  of  off-line  optimization  on  guidance  law  parameters  by  genetic  algorithm,  and  tested 
the  effect  of  guidance  parameters  under  the  conditions  of  different  initial  state  bias.  However,  from  the  simulation  results, 
the  optimization  process  costs  so  much  time  that  the  guidance  law  unfits  for  the  on-line  application.  For  more  complicated 
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situations,  a  new  integrated  guidance  and  control  design  scheme  based  on  VSC  approach  for  missile  with  terminal  impact 
angle  constraint  was  proposed  in  Wu  and  Yang  [34],  wherein  nonlinear  transformation  was  employed  to  transform  the 
mathematical  model  into  a  standard  form  suitable  for  SMC  method  design.  Meanwhile,  Gu  et  al.  [35]  presented  a  variable 
structure  guidance  law  with  constraint  on  impact  angle  that  can  be  implemented  conveniently  with  little  target  information. 
And  the  disturbances  such  as  the  estimation  errors,  the  target’s  velocity,  and  the  target’s  maneuver  acceleration  can  be  suc¬ 
cessfully  restrained  by  choosing  appropriate  guidance  parameters.  For  further  research,  a  new  scheme  of  integrated  guidance 
and  autopilot  design  was  proposed  in  Guo  and  Zhou  [36]  by  adopting  the  SMC  for  the  case  of  homing  missile  against  ground 
fixed  target  with  impact  angular  constraint.  And  the  same  problem  was  solved  using  control  method  in  Guo  and  Zhou 
[37],  Based  on  the  finite  time  convergence  stability  theory  and  the  VSC,  a  finite  time  convergent  sliding  mode  guidance 
law  with  terminal  impact  angle  constraint  was  presented  in  Zhang  et  al.  [38],  At  the  same  time,  a  3-dimensional  robust  guid¬ 
ance  law  was  developed  for  impact  angle  control  in  Hu  et  al.  [39]  by  combining  the  linear  quadratic  optimal  theory  with  the 
VSC  methodology.  Recently,  an  approach  to  impact  time  and  angle  guidance  was  presented  in  Harl  and  Balakrishnan  [40] 
through  a  combination  of  a  novel  line-of-sight  (LOS)  rate  shaping  technique  and  a  new  second-order  sliding  mode  approach. 
Then,  a  kind  of  SMC-based  guidance  law  was  proposed  in  Kumar  et  al.  [41,42]  to  intercept  stationary,  constant  velocity,  and 
maneuvering  targets  at  a  desired  impact  angle,  which  was  defined  in  terms  of  a  desired  LOS  angle.  In  addition,  an  autopilot 
was  designed  in  Kumar  et  al.  [42]  as  a  dynamic  sliding-mode  controller  that  enables  the  tracking  of  the  command  generated 
by  the  guidance  laws.  The  guidance  laws  could  intercept  the  targets  at  any  impact  angle  even  in  the  presence  of  large  initial 
heading  errors,  so  the  so-called  “all-aspect  interceptor”  could  be  achieved. 

Except  for  the  approaches  introduced  above  to  the  problem  of  impact  angle  control  for  terminal  guidance,  there  are 
many  other  studies  focusing  on  the  same  problem  based  on  different  methods.  Manchester  and  Savkin  [43]  presented  a 
precision  guidance  law  named  “Circular  Navigation  Guidance"  (CNG)  with  impact  angle  constraint  for  a  two-dimensional 
planar  intercept  based  on  the  principle  of  following  a  circular  arc  to  the  target.  As  and  different  trial,  the  impact  angle 
control  problem  was  formulated  as  an  infinite  horizon  non-linear  regulator  problem  and  solved  by  the  State-Dependent 
Riccati  Equation  (SDRE)  technique  in  Ratnoo  and  Ghose  [44,45],  Similar  study  based  on  modifying  SDRE  technique  was 
carried  out  by  Bardhan  and  Ghose  [46],  Considering  the  performance  limit  of  the  Separated  Guidance  and  Control  (SGC) 
loop,  a  concept  of  the  Integrated  Guidance  and  Control  (1GC)  was  introduced  in  Yun  and  Ryoo  [47]  and  a  new  type  of 
1GC  with  impact  angle  constraint  was  derived  using  the  Linear  Quadratic  Regulator  (LQR).  Besides,  a  cubic  spline  guid¬ 
ance  law  was  obtained  from  cubic-spline-curve-based  trajectory  using  an  inverse  method  in  Dhabale  and  Ghose  [48]  for 
intercepting  a  stationary  target  at  a  desired  impact  angle.  The  cubic  spline  trajectory  curve  expressed  the  altitude  as  a 
cubic  polynomial  of  the  downrange.  Then,  a  new  three-dimensional  guidance  law  with  high  impact  angle  constraints 
used  for  striking  ground  moving  targets  was  designed  in  Zhang  et  al.  [49]  based  on  back-stepping  control  approach. 
To  ensure  the  required  impact  angle  and  small  miss-distance,  two  step  controllers  were  designed  based  on  Lyspunov 
method.  Aiming  at  the  shortage  of  CLC  (circular  arc-straight  line-circular  arc)  path  planning,  CCC  (three  circular  arcs) 
path  planning  was  adopted  to  control  impact  time  and  impact  angle  in  Zhang  et  al.  [50],  The  latest  research  in  this  area 
can  be  found  in  Lee  et  al.  [51],  which  proposed  the  impact  angle  control  guidance  laws  with  terminal  acceleration  con¬ 
straints  for  a  stationary  or  slowly  moving  target.  These  laws  are  called  Time-to-go  Polynomial  Guidance  (TPG)  which  as¬ 
sumes  the  guidance  command  as  a  polynomial  function  of  time-to-go  and  determines  the  coefficients  of  the  guidance 
command  to  satisfy  the  specified  terminal  constraints. 

As  can  be  seen  from  the  optimal  guidance  laws,  the  guidance  problems  with  terminal  angular  constraints  are  transformed 
into  optimal  control  problems  with  terminal  constraint  conditions,  and  then  forms  the  optimal  guidance  laws  based  on  the 
OCT,  which  could  minimize  the  velocity  lost,  the  flying  time,  or  the  fuel  cost.  According  to  the  deviation  of  terminal  PNGs, 
one  or  more  terms  are  added  in  the  expression  of  the  guidance  laws  to  satisfy  the  terminal  angular  constraints.  Besides,  the 
desired  impact  angle  constraint  can  also  be  included  into  the  sliding  mode  to  form  the  sliding  mode  guidance  laws  in  more 
robust  ways.  All  of  these  guidance  laws  are  designed  not  only  to  get  a  minimum  miss-distance  but  also  to  achieve  a  desired 
impact  attitude  angle  to  acquire  better  kill  effect;  however,  the  complex  environment  the  vehicles  will  experience  and  the 
requirements  of  fuel  cost  during  the  flying  course  are  neglected.  Actually,  the  guided  weapons  will  face  up  to  harsh  condi¬ 
tions  and  typical  constraints  brought  by  dense  atmosphere  during  the  striking,  for  example,  the  heating  rate,  normal  load 
factor  and  dynamic  pressure.  However,  these  severe  conditions  have  not  been  considered  in  former  studies.  Though  the 
off-line  optimization  on  guidance  law  parameters  by  genetic  algorithm  is  studied  in  Hu  et  al.  [33],  the  whole  searching  pro¬ 
cess  casts  too  much  time  to  be  applied  on-line.  So  the  main  concern  of  this  paper  is  the  multi-constraints-based  optimization 
of  terminal  guidance  law  with  its  on-line  application. 

In  detail,  the  primary  works  and  contribution  of  this  paper  can  be  concluded  as  follows: 

1.  A  quasi-sliding  mode  guidance  law  is  design  using  boundary  layer  method  for  terminal  impact-angle-constrained 
engagement  with  stationary  target  and  the  influence  of  guidance  parameters  (i.e„  reaching  law  factor,  switching  item 
gain,  angle  error  item  factor,  and  boundary  layer  thickness)  on  the  guidance  performance  is  analyzed. 

2.  In  addition  to  the  requirement  of  impact  angle  constraint,  we  further  consider  several  hard  constraints  (i.e.,  heating  rate, 
normal  load  factor  and  dynamic  pressure)  and  control  constraints  (i.e.,  maximum  angle  of  attack  and  sideslip  angle)  dur¬ 
ing  the  striking.  Then  the  optimization  model  is  established  to  minimize  the  integrated  performance  index  consisting  of 
the  fuel  expenditure  and  velocity  lost. 
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3.  In  order  to  achieve  the  on-line  optimization  design  of  the  terminal  sliding  mode  guidance  law,  we  implement  a  fast  con¬ 
vergence  Sequential  Quadratic  Programming  (SQP)  algorithm  to  solve  the  established  optimization  problem.  Simulation 
results  show  that  the  whole  optimization  process  can  be  completed  in  a  few  seconds,  which  promotes  the  on-line  appli¬ 
cation  of  this  approach. 

The  rest  of  this  paper  is  organized  as  follows.  Section  2  establishes  the  equations  of  motion  for  the  weapon-target  engage¬ 
ment.  A  terminal  guidance  law  based  on  the  SMC  methodology  is  designed  in  Section  3  to  eliminate  the  miss-distance  and 
achieve  a  desired  impact  angle.  Then  the  influence  of  each  guidance  law  parameter  on  the  guidance  performance  is  discussed 
in  Section  4.  In  order  to  achieve  the  on-line  optimization  and  application  of  the  guidance  law,  the  optimization  model  based 
on  SQP  algorithm  is  established  in  Section  5.  Simulation  studies  are  carried  out  in  Section  6  to  investigate  the  feasibility  of 
the  SQP  algorithm  for  the  on-line  optimization  of  the  guidance  parameters.  Finally,  conclusions  and  possible  future  work  are 
discussed  in  Section  7. 

2.  Equations  of  motion 

In  order  to  simplify  the  equations  of  the  striking  situation,  the  weapon-target  engagement  model  is  divided  into  diving 
plane  and  turning  plane  as  shown  in  Fig.  1  (see  [52]).  An  earth-fixed  coordinate  system  is  defined.  The  target  is  stationary 
and  stays  at  the  origin  of  the  coordinate  system.  They-axis  is  pointed  to  the  up  and  deviates  from  the  Earth  center,  the  x-axis 
stays  in  the  horizontal  plane  and  is  pointed  to  the  weapon,  and  the  z-axis  completes  the  right-hand  system.  It  is  also  assumed 
that  the  angle  of  attack  is  small  and  its  velocity  is  constant. 

The  standard  three-dimensional  equations  of  motion  of  the  guided  weapon  over  a  flat  earth  can  be  represented  by  the 
following  nonlinear  differential  equations  (see  [52]). 

Iv  =  (Rxh  +  Pxh)/m  +  gxh, 

0  =  (Ryh  +  Pyh)/(mv)+gyh/v, 

6  =  ~(Rzh  +PZh)/(mv cos  9)  -  gzh/{ v cos  S) , 
k=  vcosOcosu, 
y  =  ysinP, 
z  =  -^cosPsin  (7, 

where  x,  y,  z  are  the  change  rate  of  position  coordinates  x,  y,  z,  respectively,  v  is  the  earth  relative  velocity  and  v  denotes  the 
acceleration.  The  flight  path  angle  and  angular  rate  are  9  and  9,  while  the  velocity  deflection  angle  and  angular  rate  are  a  and 
&.  Rxh,  Ryh  and  Rzh  are  the  aerodynamic  force  decomposition  described  in  velocity  coordinate  system  as  gxh,  gyh,  gzh  and  Pxh,  Pyh, 
PZh  represent  the  gravity  acceleration  decomposition  and  propulsion  decomposition,  respectively.  The  mass  m  is  a  constant. 

As  shown  in  Fig.  1 ,  the  motion  of  weapon  can  be  decomposed  in  diving  plane  and  turning  plane.  The  LOS  p  from  the  target 
to  the  weapon  is  defined  by  elevation  angle  XD  and  azimuth  angle  XT.  vD  and  v-r  are  the  velocity  decomposition  respectively  in 
diving  plane  and  turning  plane.  Xu  is  the  azimuth  angle  in  the  turning  plane.  The  equation  of  weapon-target  relative  move¬ 
ment  can  be  found  in  Zhao  [52] 

/  \d  =  (v/v  -  2 p/p)XD  -  pyD/p,  (2) 

\  Xtt=  (v/v  -  2p/p)Xtt  +  pyT/p, 

where  yD  and  yT  denote  the  azimuth  angle  of  velocity  in  diving  plane  and  turning  plane  respectively,  and  their  angular  rate 
are  yD  and  yT.  p  is  the  change  rate  of  LOS.  XD  and  XD  are  the  angular  rate  and  angular  acceleration  of  elevation  angle  XD,  X  n-  and 
Xjj  denote  the  angular  rate  and  angular  acceleration  of  azimuth  angle  in  the  turning  plane  Xn- 


y 


Fig.  1.  Coordinate  system  and  geometry. 
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3.  Variable  structure  guidance  law  derivation 

Due  to  the  inherent  robustness  and  relatively  simple  control  algorithm  of  the  sliding  mode  VSC  system,  the  VSC  theory  is 
gradually  applied  to  the  design  of  guidance  laws.  Most  of  the  current  researches  on  variable  structure  guidance  laws  are  con¬ 
fined  to  two-dimensional  plane,  and  how  to  select  the  guidance  parameters  also  lack  a  conclusion.  Because  the  sliding  mode 
can  be  designed  according  to  need  and  is  independent  of  parameters  and  disturbances  of  the  object,  the  desired  impact  atti¬ 
tude  can  be  included  into  the  sliding  mode  and  derive  the  variable  structure  guidance  laws  with  terminal  angular  constraint, 
which  performs  strong  robustness  on  uncertainty  factors.  Next  we  will  design  the  guidance  laws  in  diving  plane  and  turning 
plane,  respectively,  and  then  introduce  the  Quasi-Sliding  Mode  Variable  Structure  Control  (QSMVSC)  theory  to  reduce  the 
chattering  brought  by  switching. 

3.1.  Diving  plane  guidance  equation 

In  order  to  apply  the  SMC  methodology  to  the  guidance  problem,  a  switching  surface  is  to  be  chosen  first.  In  the  guidance 
problem  with  terminal  angular  constraint,  the  guidance  goal  is  transformed  into  achieving  zero  miss-distance  and  a  desired 
impact  attitude  angle,  simultaneously.  So  the  switching  surface  should  include  these  two  factors.  The  miss-distance  condi¬ 
tion  is  represented  by  LOS  rate  /.D  =  0,  and  the  terminal  impact  attitude  angular  constraint  is  satisfied  by  nulling  LOS  error 
Ad  +  }-'DF|.  So  the  switching  surface  of  this  section  is  chosen  such  that 

Si  =  Ad  +  Ai  V  ■  (Ad  +  Tdf)/ Pj  (3) 

where  the  angle  error  item  factor  Ai  is  a  positive  constant.  The  first  term  is  responsible  for  the  miss-distance  and  the  second 
term  is  responsible  for  the  impact  attitude  angle,  yDr  denotes  the  desired  impact  angle,  which  is  a  constant  since  the  target  is 
fixed. 

A  sufficient  condition  in  order  to  guarantee  the  attractiveness  of  sliding  mode  Si  =  0  is  Si  ■  Si  <  0.  The  reaching  law  is  cho¬ 
sen  as 

S,  =  -ki/Tg  ■  Si  -  £i/p  •  sgn(Si),  (4) 

where  the  reaching  law  factor  k i  and  the  switching  item  gain  £i  are  both  positive  constants.  Tg  is  the  time-to-go  and  sgn 
denotes  the  sign  function. 

Differentiating  Eq.  (3)  with  respect  to  t  and  combined  with  Eq.  (1),  the  following  holds. 

Si  =  (v/v  -  2p/p)XD  -  pyD/p  +  Ai  v(lD  +  yDF)/p  +  Ai  vXD/p  -  Ai  vp(XD  +  yDF)/p2.  (5) 

Comparing  Eq.  (4)  with  Eq.  (5),  and  considering  Tg  =  -p/p,  the  diving  plane  guidance  equation  can  be  obtained. 

7d  =  ~Tg{{v/v  +  (k,  +  2 )/Tg  +  Ai v/p)XD  +  \Xxv/p  +  Xxv(kx  +  1  )/(pTg)]  ■  (AD  +  yDF)  +  £i/p  ■  sgn(Si)}.  (6) 


3.2.  Turning  plane  guidance  equation 

In  the  turning  plane,  we  only  consider  the  condition  Xjr(tf)  =  0  at  final  time  tf.  Similarly,  the  switching  surface  is  chosen  as 
S2  =  Arr  and  the  reaching  law  is  chosen  as  S2  =  -k2/Tg  Si  -  e2/p  •  sgn(S2).  Considering  Arr  =  AT  cos  AD,  the  turning  plane 
guidance  equation  can  also  be  easily  obtained. 

yT  =  Tg{[v/v  +  (k2  +  2)/Tg]Ar  cos  AD  +  s2/p  ■  sgn(S2)}.  (7) 


3.3.  Guidance  law  with  quasi-sliding  mode  control 


The  variable  structure  guidance  law  performs  good  stability  at  the  expense  of  high  frequency  chattering,  which  can  re¬ 
duce  the  precision  of  the  control  system,  increase  the  fuel  cost  and  destroy  the  performance  of  the  system.  This  problem 
can  be  solved  by  replacing  the  sign  function  with  the  saturation  function  (see  [53]) 

M  s>  a, 

sat(S)  =  I  S/A  |:S|  <  A,  (8) 

[-1  5  <  —A, 


where  the  thickness  of  boundary  layer  A  is  a  positive  constant.  Then,  expressions  (6)  and  (7)  become  (see  [33]) 

(  h  =  -Tg{(v/v+  (ki  +  2)/Tg  +  Ai  v/p)XD  +  [A,  v/p  +  A,  p(kt  +  l)/(pTg)]  ■  (AD  +  yDF)  +  £i/p  •  sat(Si)}, 

1  yr  =  Tg{[v/v  +  {k2  +  2)/Tg]AT  cos  AD  +  s2/p  ■  saf(S2)}, 
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where  k,  X,  £  and  A  with  subscripts  1  and  2  denote  guidance  parameters  of  diving  and  turning  plane,  respectively.  This  is  the 
terminal  guidance  law  based  on  QSMVSC  methodology  to  minimize  the  miss-distance  and  satisfy  the  impact  angle  con¬ 
straint.  This  guidance  law  can  weaken  the  chattering  brought  by  switching  and  is  not  sensitive  to  external  disturbances. 

4.  Influence  of  guidance  parameters 

The  variable  structure  guidance  law  includes  four  guidance  parameters:  k,  s,  X  and  A.  It  is  necessary  to  analyze  the  influ¬ 
ence  of  the  four  parameters  on  the  guidance  performance  before  the  optimization.  It  is  known  that  the  movement  of  ground 
striking  weapon  is  mainly  in  the  diving  plane  and  the  side  guidance  has  less  impact  on  the  precision  of  miss-distance  and 
terminal  attitude  angle,  so  only  discussing  the  influence  of  ki,  X%,  £i  and  A,  is  acceptable. 

In  order  to  analyze  the  influence  of  the  four  parameters  and  investigate  the  properties  of  the  guidance  law  for  impact 
angle  control,  three-dimensional  nonlinear  simulations  with  the  initial  conditions  shown  in  Table  1  are  performed.  The  de¬ 
sired  impact  angle  is  chosen  as  yor  =  -90°. 

4.1.  Reaching  law  factor,  k, 

From  Eq.  (4),  we  can  see  that  the  factor  k,  influences  the  reaching  rate  of  the  sliding  mode.  Then  we  chose  fit,  Xu  A,  as  the 
standard  values  £t  =ss  1,  Xi  =  1,  A]  =  0.001,  and  discuss  the  effect  of  the  reaching  law  factor  k,  on  the  guidance  performance. 
Table  2  shows  the  simulation  results. 

From  Table  2,  we  can  find  that  the  terminal  velocity,  angle  of  attack,  impact  angle  error  and  the  Circular  Error  Probable 
(CEP)  change  very  small  with  different  fcis;  but  the  performance  index J,  which  represents  the  fuel  expenditure  and  velocity 
lost,  decreases  with  increasing  k,.  The  performance  index  J  is  defined  by 


(Cirf  +  C2(X2)dt, 


(10) 


where  ny  is  the  non-dimensional  normal  load  factor  in  g  and  a  is  the  angle  of  attack  in  deg.  Ci  and  C2  are  the  weighting  coef¬ 
ficients  that  are  assumed  to  be  1  and  10,  respectively. 

4.2.  Switching  item  gain,  £, 

From  Eq.  (4),  we  can  find  that  the  gain  Sj  also  affects  the  reaching  rate  of  the  sliding  mode.  Then  we  chose  Iq,  Xt ,  A\  as 
kj  =  1,  Xj  =  1,  A\  -  0.001,  and  analyze  the  impact  of  the  switching  item  gain  St  on  the  guidance  performance.  The  simulation 
results  are  as  shown  in  Table  3. 

Table  3  shows  that  the  three  trajectories  are  basically  coincident  and  the  terminal  parameters  also  behave  tiny  changes 
with  changing  £i ;  however,  with  increasing  £,,  the  performance  index  J  diminishes,  and  the  magnitude  of  normal  load  factor 
and  angle  of  attack  change  intensely. 

4.3.  Angle  error  item  factor,  Xi 

Eq.  (3)  indicates  that  the  factor  Xi  determines  the  weight  of  angle  error  item  in  switching  function  S.  Choosing  ki,  elp  Ai  as 
the  standard  values  ki  =  1,  E,  =  1,  /I  n  =  0.001,  and  discussing  the  influence  of  angle  error  item  factor  X,  on  the  guidance  per¬ 
formance.  The  simulation  results  are  reported  in  Table  4. 

It  is  evident  from  Table  4  that  with  increasing  Xi ,  the  CEP  and  terminal  angle  of  attack  change  a  little,  and  the  impact  angle 
error  and  terminal  velocity  keep  nearly  invariable.  The  numerical  value  of  performance  index  J  varies  enormously.  So  the 
factor  k,  plays  a  more  important  role  in  guidance  performance  compared  with  fc,  and  £,. 

4.4  Boundary  layer  thickness,  A  j 

According  to  the  quasi-sliding  mode  control  theory,  the  thickness  of  boundary  layer  Aj  has  an  influence  on  chattering  and 
guidance  performance.  Similarly,  we  chose  k,  and  Xi  as  the  standard  values,  and  £i  =  10.  Then  analyze  the  impact  of  the 
boundary  layer  thickness  Ai  on  the  guidance  performance.  The  simulation  results  are  as  follows. 


Initial  position  (x0,yo,z0) 
Initial  velocity  ( va ) 


(10,16,-10)  km 
340  m/s 


Initial  heading  angle  (e0) 
Initial  deflection  angle(co) 


-120° 
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Table  2 

Comparisons  of  the  effects  for  different  fcis  ()>df  =  -90°). 


fcl 

CEP(  rn) 

0/(°) 

VfW  s) 

J 

0.1 

0.0798 

-90.1565 

-0.1750 

677.2032 

749.3447 

1.0 

0.2097 

-89.9250 

-0.1725 

681.0937 

746.5281 

10.0 

0.6126 

-89.2406 

-0.0388 

680.5330 

737.7190 

Table  3 

Comparisons  of 

the  effects  for  different  s,  s  (yDF  = 

■  -90"). 

E, 

CEP  (m) 

0/(°) 

%(") 

vf  (m/s) 

J 

0.1 

0.0441 

-90.1663 

-0.1769 

677.2492 

750.3775 

1.0 

0.2097 

-89.9250 

-0.1725 

681.0937 

746.5281 

10.0 

0.1883 

-89.9750 

-0.1791 

680.3869 

737.9296 

Table  4 

Comparisons  of  the  effects  for  different  A,s  (7df  = 

=  -90"). 

M 

CEP  (m) 

%(“) 

«/(-) 

ijr(m/s) 

j 

0.1 

0.9037 

-89.4823 

-0.2938 

676.5685 

707.2557 

1.0 

0.2097 

-89.9250 

-0.1725 

681.0937 

746.5281 

10.0 

0.0446 

-90.0260 

-0.0081 

681.1353 

1201.5332 

Table  5 

Comparisons  of  the  effects  for  different  A\ s  (yDF 

=  -90"). 

2fi 

CEP  (m) 

S/(°) 

®r« 

tytm/s) 

J 

0.0001 

0.8566 

-89.5853 

-0.1103 

680.0898 

762.3297 

0.001 

0.1883 

-89.9750 

-0.1791 

680.3869 

737.9296 

0.01 

0.3083 

-90.2705 

-0.1864 

681.0053 

745.1327 

Table  5  shows  that,  the  guiding  precision  and  controlling  performance  change  very  small  with  different  z^s,  but  under¬ 
size  zh  can  cause  severe  chattering,  which  can  reduce  the  precision  of  the  guidance  system. 

In  general,  the  selection  of  guidance  parameters  should  guarantee  the  high  guiding  precision  first.  To  achieve  better  strik¬ 
ing  effect  and  avoid  the  saturated  command,  minimizing  the  miss-distance  and  terminal  normal  load  factor  is  required.  In 
addition,  the  fuel  cost  and  the  constraints  of  heating  rate,  normal  load  factor  and  dynamic  pressure  should  also  be  consid¬ 
ered.  So  the  parameters  of  variable  structure  guidance  law  should  be  designed  and  optimized  under  above  considerations. 


5.  On-line  optimization  of  guidance  law  with  SQP  algorithm 

The  feasibility  of  off-line  optimization  on  guidance  law  parameters  by  genetic  algorithm  is  discussed  in  Hu  et  al.  [33];  but 
it  costs  so  much  time  and  limits  the  on-board  application  of  the  guidance  law.  In  order  to  accomplish  the  goal  of  on-line 
application  for  the  ground  striking  weapon,  it  is  necessary  for  us  to  redesign  the  optimal  guidance  law  parameters  under 
multiple  constraints  with  other  efficient  optimization  methods. 

Many  optimization  problems  in  aeronautics  and  astronautics  can  be  formulated  as  optimal  control  problems.  Stryk  and 
Bulirsch  [54]  gives  a  brief  list  of  commonly  used  direct  and  indirect  efficient  methods  for  the  numerical  solution  of  these 
problems.  For  some  special  weakly  nonlinear  low  dimensional  systems  the  solution  can  be  derived  analytically  from  the  suf¬ 
ficient  and  necessary  conditions  of  optimality;  but  to  obtain  a  solution  of  dynamic  systems  described  by  strongly  nonlinear 
differential  equations,  it  is  necessary  to  use  numerical  methods. 

Stochastic  approaches  such  as  Genetic  Algorithm  (GA),  Particle  Swarm  Optimization  (PSO)  and  Simulated  Annealing  (SA) 
algorithms  are  used  for  many  kinds  of  optimization  problems  (see  [55-58]).  Many  other  methods  are  also  proposed  to  im¬ 
prove  the  optimization  performance  (see  [59,60,56,57]).  In  recent  years,  general-purpose  SQP  methods  have  been  consider¬ 
ably  developed  and  applied  to  trajectory  optimization  and  other  optimization  problems  with  linear  and  nonlinear 
constraints  (see  [61-64]).  Some  practical  and  theoretical  aspects  of  applying  SQP  methods  to  optimal  control  problems 
are  discussed  in  [65].  These  methods  can  be  shown  to  converge  to  a  solution  under  very  mild  conditions  on  the  problem. 
Terminal  guidance  law  on-line  optimization  problem  is  a  typical  nonlinearly  constrained  optimization  problem,  which 
can  be  solved  by  SQP  methods. 

The  performance  index  is  defined  by  Eq.  (10)  which  represents  the  fuel  expenditure  and  the  velocity  lost  during  the  strik¬ 
ing  course.  At  the  same  time,  the  guidance  law  parameters  should  also  be  selected  to  satisfy  the  typical  constraints  such  as 
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heating  rate,  normal  load  factor  and  dynamic  pressure  constraints,  and  maximum  angle  of  attack  and  sideslip  angle  con¬ 
straints.  So  our  goal  is  to  search  for  adaptable  guidance  law  parameters  which  not  only  satisfy  the  specified  constraints 
but  also  minimize  the  performance  index. 

Here  we  divide  the  total  time  interval  [0,t/]  into  n  -  1  equal  subintervals  of  length  h  each  by  the  points  {ti,  t2,  ■  ■ .  t„}  with 
^  =  0  and  tn  =  tf.  Then  minimizing  the  performance  index  is  equivalent  to  minimizing 

r=lY^(C,n2yi  +  C2xt)-h.  (11) 

When  the  step  size  h  is  fixed,  it  is  appropriate  to  minimize  the  numerical  value  of  (C,  nj.  +  C2oif)  relates  to  the  real-time  guid¬ 
ance  law  parameters  in  each  interval  instead  of  (11). 

In  each  interval  [t;,tj+i],  i  =  1,2 . n  —  1,  when  given  the  vector  of  state  variables  X,  =  (i>it  8t,  at,  x(,  yf,  z,),  the  current 

guidance  law  can  be  derived. 

j  7d  =  -Tg[{Vi-i/Vi  +  (kj  +  2)/Tg  +  hVi/p)\D  +  [Mi)i..t/p  +  Ai Vi(kt  + 1  )/(pTg)\  •  (AD  +  yDF)  +  6i/p  ■  sat(SM)}, 

[  Tt  =  Tg[[Vi-tlVi  +  (k2  +  2)/Tg]ATcos/lD+82//>  ■  sot(S2i)}, 

(12) 

where  V\  ,  is  the  velocity  rate  of  last  step,  and 
f  AD  =  arctan(yf /^xj+zj), 

[  At  =  arctan(-Zi/Xi), 

{Ad  =  Vjp, 

At  =  —v^l(pcoslp), 

P  =  y/xj+yf+zl 
Tg  =  p/vt  =  -p/p, 


{Xi  =  Vi  COS  (/  COS  (Ti, 

y,  =  Vi  sin  (),, 

Zi  =  -  Vi  cos  dj  sin  di, 


[cos  Ad  cos  At  sin  Ad  -  cos  AD  sin  lT 
-sin  Ad  cos  At  cosAd  sin  AD  sin  Ar 

sin  At  0  cosAt 

where  v(,  vn  and  Vr  are  the  velocity  decomposition  described  in  the  LOS  Cartesian  coordinate  system.  So  denotes  the  trans¬ 
form  matrix  between  the  target  coordinate  system  and  the  LOS  coordinate  system.  Then  we  obtain 

9  =  -y d/cos(At  -  a/) , 

<7=  [yr- I'D  tan( At  -  at)  sinAD]/cosAD. 

From  (1 ),  we  get  the  lift  and  side  force 

Yt  =  mvi  {8-  p/(r3Vi )  •  [xf  sin  0,  cos  a,  -  (yf  +  R0)  cos  Ot  -  z,  sin  (/  sin  (Ti]}  -  Psina^i, 

Zi  =  -mvi  cos  dr  {a  +  p/(r3Vi  cos  0,-)  ■  [x<  sin  c,  +  z<  cos  at]}  +  P  cos  an  sin  [/  , , 

where  r  is  the  radial  distance  from  the  center  of  the  earth  to  the  weapon  and  p  is  the  gravitational  constant  while  R0  denotes 
the  mean  radius  of  the  Earth.  af_i  and  ,  are  the  angle  of  attack  and  sideslip  angle  of  last  step.  The  propulsion  P  is  a 
constant. 

It  is  assumed  that  the  angle  of  attack  and  the  sideslip  angle  are  small.  Then  we  can  approximately  obtain 

«i  —  Yi/(Cy  ■  q-S), 

Pi  =  Zi/(C^q-S), 


(13) 
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Fig.  2.  The  circuit  of  on-line  optimization. 


where  q  is  the  dynamic  pressure.  S  denotes  the  reference  area.  C“  and  Cf  are  the  aerodynamic  lift  and  side  force  coefficients, 
respectively.  The  absolute  magnitude  of  angle  of  attack  and  sideslip  angle  should  be  limited  during  the  whole  striking  course 
by  their  maximum  values  amax  and  /3max. 
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Velocity  versus  time 


t(s) 


Fig.  4.  Velocity  profile  for  the  ground  striking  (yDF=  -45°). 


The  drag  is  approximated  by 
Xi  =  CxqS, 

where  CX  =  CX o  +  C“( a?  +  /??).  Cx0  and  Cx  are  aerodynamic  drag  coefficients. 

At  this  moment,  the  velocity  rate  is 

Vi  =  (P  cos  Of  cos  Pi  -  Xj)  /m  -  ju/r3  ■  [x,  cos  6,-  cos  +  (yt  +  R0)  sin  0t  -  z,-  cos  0,-  sin  er,]. 

Then  the  heating  rate  constraint,  normal  load  factor  constraint  and  dynamic  pressure  constraint  can  be  expressed  analyt¬ 
ically  by  expressions  of  the  guidance  law  parameters,  respectively. 


Qi  =  y/Pi(Vi- 1  +  Vi  ■  hf  ~f(kl ,  8j  ,  At)  <  Qmax, 

riyi  =  JYf+Zf  ■  cos  Oj  +X,  ■  sin  a,-  =f(k,  ,Si,M)<nymax, 
qi=0.5pi( Vi  ,  +  Vi-hf  =f(k,,e,,,M)  < qmax , 


(14) 


Flight  path  angle  versus  time 


t(s) 


Fig.  5.  Flight  path  angle  profile  for  the  ground  striking  (yDF  =  -45°). 
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Angle-of-attact  versus  time 


t(s) 

Fig.  6.  Angle  of  attack  profile  for  the  ground  striking  (yDF  =  -45°). 


where  Vi  ,  is  the  velocity  of  last  step.  Qmax,  nymax  and  qmax  are  the  maximum  allowable  values  of  heating  rate,  normal  load 
factor  and  dynamic  pressure,  respectively.  Besides,  the  design  variables,  i.e.  the  real-time  guidance  law  parameters,  should 
lie  in  their  upper  and  lower  bounds. 

0.01  <  k,  <  2,  0.01  <  e,  <  2,  0.1  <  Aj  <  2. 

Here  we  only  consider  three  parameters  and  choose  zli  as  the  standard  constant  because  of  its  negligible  impact  on  the  guid¬ 
ance  performance. 

Based  on  above  procession,  our  crucial  target  is  to  find  out  a  series  of  optimal  guidance  law  parameters  which  not  only 
satisfy  the  specified  inequality  constraints  as  shown  in  (13)  and  (14),  but  also  minimize  the  performance  index  in  every  step. 
In  addition,  the  design  variables  should  not  exceed  the  upper  and  lower  bounds. 

The  flow  chart  of  on-line  optimization  is  shown  in  Fig.  2. 


Normal  load  versus  time 


t(s) 


Fig.  7.  Normal  load  factor  profile  for  the  ground  striking  (yDF=  -45°). 
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Reaching  law  factor  versus  time 


t(s) 

Fig.  8.  The  optimal  history  of  It,  for  the  guidance  law  (yDF  =  -45°). 


6.  Simulation  studies 

The  desired  impact  angle  is  chosen  as  yDF  =  -45°  and  yDF  =  -90°,  respectively,  to  investigate  the  properties  of  the  SQP 
method  for  the  guidance  law  optimization.  The  maximum  allowable  values  of  these  inequality  constraints  are  set  at 
<*max  =  6°,  Anax  =  6°,  Q,max  =  1 , 400, 000  kW/m2,  Hymax  =  12  g  and  qmax  =  350,000  N/m2.  The  initial  conditions  are  given  in  Ta¬ 
ble  1.  This  simulation  task  is  carried  on  a  PC  with  2G  memory  and  Intel(R)  Core(TM)  2  Duo  CPU.  The  version  of  MATLAB  is 
R2012a. 


6.1.  Desired  impact  angle  yor  =  -45° 

The  simulation  results  are  presented  by  Figs.  3-10  and  Table  6. 


Switching  item  gain  versus  time 


t(s) 


Fig.  9.  The  optimal  history  of  n  for  the  guidance  law  (yDF  =  -45°). 
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Angle  error  item  factor  versus  time 


- -uy-A/r- 

lJIP 

iy 

0  5  10  15  20  25 


t(s) 


Fig.  10.  The  optimal  history  of  7,  for  the  guidance  law  ( yDP  =  -45°). 


Table  6 

Simulation  results  for  ground  striking  (yDF=  -45°). 

CEP  (m) Bf(°) Of« (m/s) ] 

1.5639  -44.8619  0.2287  651.7479  671.9715 


6.2.  Desired  impact  angle  yDF  =  —90° 

The  simulation  results  are  presented  by  Figs.  11-18  and  Table  7. 

Figs.  3  and  11  present  the  three-dimensional  striking  trajectories  and  their  ground  tracks.  Figs.  4  and  12  show  the  time- 
histories  of  velocity.  Figs.  5  and  13  depict  the  flight  path  angle  curves  from  the  initial  value  to  the  desired  impact  angle.  Figs. 


Three-dimensional  trajectory 


Fig.  11.  Three-dimensional  trajectory  for  the  ground  striking  (yDF  =  -90°). 
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6  and  14  represent  the  angle  of  attack  histories  and  the  corresponding  normal  load  factor  profiles  are  shown  by  Figs.  7  and 
15,  respectively.  By  observing  Fig.  14  with  a  desired  impact  angle  of  -90°,  the  values  of  angle  of  attack  remain  negative 
through  the  whole  striking  course  to  make  the  trajectory  more  precipitous,  so  that  the  flight  path  angle  can  be  reduced 
quickly.  However,  in  the  case  of  —45°  desired  impact  angle,  as  can  be  seen  from  Fig.  6,  the  magnitude  of  angle  of  attack  main¬ 
tains  negative  duding  most  part  of  the  flying  course  and  then  turns  to  positive  to  meet  the  requirement  of  the  desired  impact 
angle.  Besides,  with  the  selecting  of  the  optimal  guidance  law  parameters,  small  daps  of  the  angle  of  attack  emerge  during 
the  initial  stage  of  the  striking,  which  can  also  be  observed  from  the  history  of  normal  load  factor  shown  in  Figs.  7  and  15. 
Figs.  8-10  and  16-18  show  the  optimal  profiles  of  guidance  parameters  based  on  the  SQP  algorithm. 

Simulation  results  for  ground  striking  are  given  in  Tables  6  and  7.  By  these  results,  we  can  find  that  the  SQP  method  for 
on-line  optimization  design  of  guidance  law  with  multiple  constraints  is  well  estimated.  Furthermore,  when  the  step  size  is 
chosen  as  0.5  s,  the  whole  optimization  process  only  costs  less  than  5  s,  while  the  off-line  optimization  by  genetic  algorithm 
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Angle-of-attact  versus  time 


t(s) 

Fig.  14.  Angle  of  attack  profile  for  the  ground  striking  (yDF=  -90°). 


Normal  load  versus  time 


t(s) 

Fig.  15.  Normal  load  factor  profile  for  the  ground  striking  ()>df=  -90°). 


discussed  in  Hu  et  al.  [33]  cost  about  20  min.  So  the  optimization  strategy  proposed  in  this  paper  satisfies  the  requirement  of 
real-time  optimization  and  on-line  application. 

6.3.  Simulation  results  for  different  impact  angles 

An  optimization  procedure  of  SMC-based  terminal  guidance  law  for  impact  angle  constrained  engagement  with  station¬ 
ary  target  is  proposed  and  two  cases  are  studies  to  illustrate  the  feasibility  of  this  method.  However,  it  is  still  not  enough  to 
demonstrate  and  address  the  performance  of  this  new  approach.  In  this  subsection,  the  magnitude  of  impact  angle  is  as¬ 
sumed  to  range  from  0  to  -90°,  and  the  corresponding  miss-distances,  impact  angles,  terminal  velocities  and  angles  of  attack 
are  reported  in  Table  8. 
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Reaching  law  factor  versus  time 


t(s) 


Fig.  16.  The  optimal  history  of  it,  for  the  guidance  law  (yDF=  -90°). 


Switching  item  gain  versus  time 


I'H 

1 

. : . i . IV 

n  I  : 

0  5  10  15  20  25 


t(s) 

Fig.  17.  The  optimal  history  of  e,  for  the  guidance  law  (yDP  =  -90°). 


As  can  be  seen  from  Table  8,  the  air-to-ground  weapon  can  achieve  nearly  every  desired  impact  angle  ranging  from  0  to 
-90°  with  high  performance  after  using  the  sliding  mode  terminal  guidance  law  and  on-line  optimization  procedure.  Both 
the  miss-distances  and  the  terminal  angles  of  attack  are  small,  and  the  values  of  terminal  heading  angle  are  almost  the  same 
as  each  desired  impact  angle.  Furthermore,  the  optimal  solutions  of  guidance  parameters  under  each  case  can  be  obtained 
within  several  seconds  based  on  SQP  method. 


6.4.  Simulation  results  under  initial  errors 

The  air-to-ground  guided  weapons  are  usually  carried  by  airplane  or  even  larger  transporter  to  specified  position  and 
launched  to  carry  out  the  striking.  However,  the  accurate  toss  condition  is  difficult  to  reach  because  of  unknown  distur- 
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Angle  error  item  factor  versus  time 
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Fig.  18.  The  optim 

Table  7 

Simulation  results  for  ground  striking  (yDF  =  -90°). 

10 

t(s) 

il  history  of  7i  for  the  j 

15 

guidance  law  (yDj 

20 

f  =  —90°). 

25 

CEP  (m) 

0/(°) 

Vf  (m/s) 

j 

0.0718 

-90.1115 

-0.0912 

679.1426 

664.5976 

Table  8 

Simulation  results  for  different  desired  impact  angles. 

•in f  (°) 

CEP  (m) 

0ff»i 

«/(“) 

F/(m/s) 

0 

2.8187 

-0.6613 

0.6185 

603.0048 

-15 

2.5522 

-14.5442 

0.5148 

618.9391 

-30 

2.1803 

-29.6635 

0.3802 

635.8122 

-45 

1.5639 

-44.8619 

0.2287 

651.7479 

-60 

1.1970 

-59.9764 

0.0722 

664.8032 

-75 

1.0399 

-75.0553 

-0.0701 

674.2781 

-80 

0.6916 

-80.0542 

-0.1022 

675.5686 

-85 

0.5911 

-85.1164 

-0.1513 

677.7938 

-90 

0.0718 

-90.1115 

-0.0912 

679.1426 

Table  9 

Simulation  results  under  initial 

deviations  (yDr  =  -90°). 

Initial  deviations 

CEP  (m) 

ffj 

«f(’> 

if  (m/s) 

Initial  height  h0 

+5% 

0.1628 

-90.1144 

-0.1580 

677.9521 

-5% 

0.1940 

-90.1393 

-0.1894 

676.4697 

Launching  range  p0 

+5% 

0.1834 

-90.1234 

-0.1769 

677.3472 

-5% 

0.2065 

-89.9281 

-0.1688 

681.0739 

Initial  velocity  ifo 

+5% 

0.1789 

-90.1237 

-0.1733 

677.3064 

-5% 

0.2176 

-89.9360 

-0.1718 

680.8195 

Initial  heading  angle  fl0 

-5° 

0.1928 

-90.1193 

-0.1837 

677.2561 

Initial  deflection  angle  <t0 

+5° 

0.0591 

-90.1161 

-0.1898 

679.3136 

-5° 

0.0589 

-90.1157 

-0.1901 

679.3184 

bances  such  as  navigation  errors,  atmospheric  disturbances  and  wind  interference.  A  guidance  law  with  strong  robustness 
should  be  insensitive  to  these  uncertain  factors.  In  order  to  validate  the  robustness  of  the  terminal  guidance  law  and  its  opti- 
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mization  approach,  the  effect  of  initial  condition  deviation  on  guidance  performance  is  studied  and  the  simulation  results  are 
reported  in  Table  9.  The  desired  impact  angle  is  chosen  as  yDF  =  -90°. 

According  to  the  numerical  simulation  results,  the  proposed  optimization  design  of  the  SMC-based  guidance  law  for 
stationary  target  offers  good  robustness  and  becomes  insensitive  to  uncertainties.  The  miss-distances,  impact  angles  and 
terminal  angles  of  attack  under  different  initial  deviations  approximate  to  the  results  of  ideal  condition  without  any 
disturbances. 

7.  Conclusions 

According  to  the  fact  that  many  kinds  of  air-to-ground  guided  weapons  are  hoped  not  only  to  get  a  minimum  miss-dis¬ 
tance  but  also  to  achieve  a  desired  impact  angle  to  acquire  better  kill  effect,  a  terminal  guidance  law  based  on  the  SMC 
methodology  was  designed  to  achieve  this  goal.  This  paper  studies  the  on-line  optimization  of  the  terminal  guidance  law 
using  SQP  algorithm  and  obtains  the  optimal  history  of  the  guidance  parameters  corresponding  to  specified  constraints 
and  performance  index. 

First,  the  weapon-target  engagement  model  was  divided  into  two  separate  planes  (i.e.,  the  diving  plane  and  the  turning 
plane)  and  the  equations  of  the  strildng  situation  were  derived.  Then  a  terminal  guidance  law  was  designed  using  SMC 
theory  to  null  the  miss-distance  and  achieve  a  desired  impact  angle,  simultaneously.  The  chattering  of  SMC  guidance  was 
reduced  by  the  boundary  layer  approach  and  formed  the  Quasi-SMC  guidance  law.  The  influence  of  the  four  guidance  param¬ 
eters  on  the  guidance  performance  was  discussed  as  well.  Based  on  the  consideration  of  the  real-time  requirement  and  the 
complex  environment  of  condensed  atmosphere,  the  optimization  model  of  the  guidance  law  was  established  to  satisfy  spec¬ 
ified  constraints  (i.e.,  heating  rate,  normal  load  factor,  and  dynamic  pressure)  and  minimize  the  performance  index,  which 
represents  the  fuel  expenditure  and  velocity  lost.  Finally,  this  optimization  model  was  solved  by  SQP  algorithm  and  its  fea¬ 
sibility  was  demonstrated  through  simulation  studies  under  different  scenarios. 

The  new  approach  discussed  in  this  paper  serves  as  a  useful  and  efficient  mean  on  the  design  and  optimization  of  the 
guidance  system  for  the  anti-tank  or  anti-ship  missiles,  the  precision  guided  penetrating  bombs  and  some  other  guided 
weapons.  According  to  the  simulation  results,  the  proposed  on-line  optimization  design  of  the  sliding  mode  guidance  law 
not  only  satisfies  the  specified  constraints,  but  also  minimizes  the  fuel  cost  during  the  whole  striking  course.  What  is  more 
important  is  the  optimization  process  can  be  completed  in  a  few  seconds,  which  is  suitable  for  on-board  application  of 
ground  striking.  Future  work  in  this  area  could  focus  on  developing  new  optimization  strategies  for  the  design  of  guidance 
laws.  The  proposed  new  approach  could  be  extended  to  solve  the  optimization  problems  of  new  guidance  laws  under 
different  engagement  scenarios,  for  example,  the  high  order  sliding  mode  guidance  laws  for  air-to-ground  guided  weapons, 
the  adaptive  sliding  mode  guidance  laws  for  air-to-air  missiles,  and  the  fuzzy  sliding  mode  guidance  laws  for  anti-ship  mis¬ 
siles,  and  so  on. 
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